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Abstract. We have discovered molecular hydrogen in a 
fourth quasar damped Lya system (hereafter DLA). The 
UVES spectrograph on the 8.2m ESO Kueyen telescope 
has allowed the detection of H2 in gas with low metallicity, 
Z/Zq ~ 10~^, and high neutral hydrogen column density, 
7V(Hi) ~ 2.6 X 10^^ cm~2, at redshift Zabs = 3.3901 toward 
QSO 0000-2620. The measured H2 fractional abundance 
of /(H2) ~ 4 X 10~® is lower than a typical value for 
Galactic interstellar clouds of high iV(Hi) column density 
by a factor of (2 — 3) x 10^. Since H2 molecules are formed 
efficiently on dust grains, it implies that the dust conden- 
sation in this DLA is negligible, and hence the abundances 
derived from metal absorption lines are the actual ones. 
The obtained /(H2) value leads to an estimate of the dust 
number density of (nd)DLA ^ lO^'^ (?T.d)isM, which is con- 
sistent with the dust-to-gas ratio k < 1.6 x 10^'^ derived 
independently from the [Cr/Zn] and [Fe/Zn] ratios. 

Key words: line: profiles - ISM: molecules - quasars: 
absorption lines - quasars: individual: QOOOO-2620 



1. Introduction 

In this work we continue our chemical composition anal- 
ysis of the Zabs = 3.3901 damped Lya system from the 
spectrum of the bright quasar QOOOO-2620 (V = 17.5, 
Zem = 4.108) observed with the UVES spectrograph 
(Dekker et al. 2000) at the 8.2m ESO Kueyen telescope 
during the first commissioning in October 1999 (Molaro 
et al. 2000, hereafter Paper I). 

Paper I was mainly concerned with precise mea- 
surements of the oxygen and zinc abundances. The ob- 
tained O and Zn metallicities together with other a- 
chain and iron-peak element abundances present a pat- 
tern of a galactic chemical composition different from that 
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of the Milky Way. In particular, it was found that the 
relative abundances of a-chain and iron-peak elements, 
[0,Si,S/Cr,Fe,Zn], are significantly lower than analogous 
ratios in Galactic stars with comparable metallicities. The 
measured Zn abundance, [Zn/H] — — 2.07±0.1 dexQ, is the 
lowest among DLAs, which implies that the Zabs — 3.3901 
intervening galaxy is in the early stages of its chemical 
evolution. This suggestion is consistent with the similar 
abundances found for the volatile (Zn) and refractory (Cr 
and Fe) elements, [Cr/H] = -1.99 ±0.1 dex and [Fe/H] = 
—2.04 ±0.1 dex, which means, in turn, that the Cr and Fe 
atoms may not be locked up in dust grains and that the 
dust-to-gas ratio in this DLA may be much lower if com- 
pared with Galactic diffuse clouds with the same neutral 
hydrogen column densities, A^(IIi) ~ 2.5 x 10^^ cm^^. 

The low dust-to-gas ratio in this system suggested in 
Paper I is consistent with the absence of significant H2 
absorption at Zabs — 3.3901 in the Lyman B^E+ <— X^I]+ 
and Werner C^n± ^ X^T.+ bands (Levshakov et al. 1992, 
hereafter LCFB) as well as with the absence of CO lines 
from the fourth positive band A^II ^ X^I]+ (Lu et al. 
1999). 

It should be noted that in the Galactic interstellar 
medium, sightlines with 7V(Hi) > 10^^ cm~^ show high 
fractional abundances of hydrogen in molecular form, 
/(H2) - 0.1 (Savage et al. 1977), where /(H2) = 
2iV(H2)/iV(H) with A^(H2) and A^(H) being the total col- 
umn densities of H2 and H, respectively. Such molecu- 
lar clouds with A^(H2) > 10"'^^ cm~^ provide an effective 
self-shielding in the Lyman and Werner bands preventing 
UV photons from penetrating the interior of the clouds 
and thus lowering photodissociation rates of H2 and CO 
molecules. 

In contrast with the ISM diffuse clouds, molecular hy- 
drogen in DLAs shows different behavior. Namely, we do 
not observe a pronounced transition from low to high 
molecular fractions at total hydrogen column density near 
5 X 10^° cm~^ as found in the Milky Way disk by Savage et 
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al. (1977). The first H2 system identified at Zabs = 2.8112 
toward Q0528-250 by Levshakov & Varshalovich (1985) 
gives very low amount of H2, /(H2) c± 5.4 x 10~^, while 
7V(H i) ~ 2.2 X 10^^ cm^^ (see also Srianand & Petitjean 
1998 and references cited therein). However, the second 
and the third H2 systems found at z^bs = 1-97 toward 
Q0013-004 (Ge & Bechtold 1997) and at Zabs = 2.34 to- 
ward Q1232+0815 (Ge et al. 2000) show high H2 abun- 
dances, /(H2) 0.22 and 0.07, respectively (the corre- 
sponding neutral hydrogen column densities are equal to 
6.4 X 10^° cm-2 and 8.0 x lO^o cm^^). Other DLAs have 
led so far only to upper limits, /(H2) ^ 2 x 10"^, with the 
lowest bound of 9 x 10"'^ at Za.bs = 2.47 toward Q1223+17 
(Bechtold 1999). 

LCFB derived an upper limit of /(H2) < 3 x 10"*^ 
for the DLA at Zabs = 3.3901 toward QOOOO 2620. The 
correlation between the H2 fractional abundance and the 
[Cr/Zn] ratio found by Ge & Bechtold (1999; see their 
Fig. 1) suggests that /(H2) < lO"'^ if [Cr/Zn] - 0, as 
measured in Paper I. This motivated us to search for the 
H2 transitions in the UVES spectrum already discussed in 
Paper I. We show below that the S/N of this spectrum is 
indeed adequate to measure H2, thus making this absorber 
the fourth one for which such a measure is available. 



2. Data analysis and Results 

2.1. Observations and the H2 identification 

Spectroscopic observations of QOOOO-2620 are described 

in detail in Paper I. The spectrum was obtained with the 
rms uncertainty of the wavelength calibration 6X < 0.6 
km s~^, the velocity resolution of FWHM ~ 6 km s~^ 
(the corresponding bin size equals 2.4 km s"'^) and the 
signal-to-noise ratio of S/N ~ 40 (per pixel) in the range 
AA = 4605 - 4615 A which allows us to detect the L(4- 
O)R(l) line at the expected position, Aobs = 4609.4 A. 
The absence of this line in the 1 A resolution spectrum 
obtained with the Multiple Mirror Telescope (MMT) in 
1988 (LCFB) has led to the limit W.cst < 114 mA (3cr) 
for the equivalent width of the L(4-0)R(1) fine. The UVES 
spectrum shown in Fig. la provides Wrest — 6 mA. The 
local continuum level in the vicinity of the L(4-0)R(1) line 
was additionally controlled by the optical depth value of 
the Lyman-limit discontinuity (tc = 0.003) from the metal 
system at Zabs = 4.061 (Savaglio et al. 1996). This is il- 
lustrated in Fig. 2 where the estimated continuum level 
between 4555 A and 4660 A is shown by the dotted line. 
The vertical line at A = 4614.39 A marks the Lyman- 
limit position for the ^abs = 4.061 system. The expected 
optical depth in the continuum at A = 4609 A is equal to 
T4609 = Tc (4609/4614)3 ~ 0.003. 

We have thoroughly investigated all possible coinci- 
dences of the L(4-0)R(1) position with metal absorptions 
from other systems observed toward QOOOO-2620. Only 
one such case was found. The position of the C iiA903.6235 



_> 

Pi 



L(4-0) R(l) 




11 til,' 
1 1 - 



Nj=i = 3.9 X 10'^ cm"^ b = 8.8 km/s 
Nj^i = 3.0 X 10^^ cm^ b = 8.0 km/s 

Nj^i =4.8x10^^ cm'^ b= 11.0 km/s 

1 I 1 1 1 I 1 1 1 I 1 1 1 I ■_ 



-40 



-20 







20 



40 



Relative Velocity (km / s) 




Nj=i/10" (cm"^) 

Fig. 1. (a) - Observed normalized intensities (dots with 

error bars) vs radial velocities given relative to Zabs = 
3.390127. The rest frame wavelength of the L(4-0)R(1) 
line and its oscillator strength (1049.958 A and 0.0160, re- 
spectively) were taken from Morton & Dinerstein (1976). 
Smooth lines are the synthetic spectra obtained from the 
minimization as described in the text, (b) - Confidence 
region in the 'Afj=i(H2) — b' plane. The contour above the 
dashed horizontal line represents 70% accepted confidence 
level. The horizontal line restricts 6- values which are less 
than the lowest Icr boundary for the 6-value of the heav- 
iest element from the Zabs = 3.3901 system, 6znii > 8.0 
km s~^ (Paper I). The cross marks the point of maximum 
likelihood for the best fit shown by the solid line in panel 



(/abs = 0.168) line from the metal system at Zabs = 
4.10106 detected by SavagUo et al. (1997) differs from 
the center of the L(4-0)R(1) line by only —0.7 km s~^. 
If the feature observed (Wobs = 26 mA) were due to 
this C II absorption, we would obtain W^est = 5 mA and 
iV(Cil) = 4.27 X 10^^ cm~^ using the Hnear part of the 
curve of growth. This would imply an abundance [C/H] 
~ -1-1.68 which is unlikely, due to the high redshift of the 
absorber. We remark that even assuming a solar carbon 
abundance, [C/H] = 0, and the unrealistic suggestion that 
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Fig. 2. A portion of the Q 0000 2620 spectrum in the vicinity of the H2 L(4-0)R(1) Hne (Aobs = 4609.4 A). The 
Lyman-Umit discontinuity of the metal system at ^;abs = 4.061 (tc = 0.003) is marked by the vertical line (see text for 
details) 



all carbon is in the C II form (Savaglio et al. have detected 
Civ in this system which is not a DLA) we find that 
Wobs(C 11) = 0.5 mA. This result implies that the C 11 ab- 
sorption at 2abs = 4.10106 does not affect significantly the 
H2 profile. 

The identification of other H2 transitions from the Ly- 
man and Werner bands is hampered by blending with nu- 
merous Lya forest lines. Therefore, to obtain quantitative 
estimates of the total molecular hydrogen column den- 
sity and the excitation temperature, Tex, we choose small 
continuum 'windows' which are less obscured by the Lya 
forest absorption features and where some H2 lines are 
expected. The following is a description of the adopted 
procedure. 

2.2. Column densities and the excitation temperature 

The inferred H2 column density, Niji^) = ^j^Ji de- 
pends on the populations of the first rotational levels 
(J = 0, 1, 2, ...) in the ground state, Tex, and the Doppler 
6-parameter. In thermal equilibrium, the ratio of the col- 
umn densities Nj is given by 



No 



gj_ 

.90 



■ exp 



B„J(J + 1) 



(1) 



where the statistical weight of a level J is 3(2 J-|- 1) for odd 
J and (2 J-l- 1) for even J (the ortho- and para-H2, respec- 
tively), and the constant By = 85.36 K for the vibrational 
ground state. 

The central optical depth, tq, for the resolved absorp- 
tion line is determined by (e.g., Spitzer 1978) 



TO = 1.497 X 10-15 Ao/absX ' 





(2) 



where /abs is the oscillator strength of the transition, Aq 
its wavelength in A, N the total column density in cm"^, 
and b the Doppler parameter in km s-^. For the H2 lines, 
wavelengths and oscillator strengths were taken from the 
list of Morton & Dincrstein (1976). 

To estimate the ortho-H2 column density, Ni , and the 
6-parameter, we fit theoretical Voigt profiles to the ob- 
served intensities via minimization using a fixed value 
of Zabs = 3.390127 which corresponds to the redshift of 
the isolated NiiiA1709 absorption line (Paper I). During 
the fitting procedure the theoretical profiles arc convolved 
with a Gaussian of the instrumental width of 6 km s^^. 
We tried to find a satisfactory result with the reduced 

per degree of freedom of xlun < X^,o.30 (where xl,o.3o 
is the expected value for v degrees of freedom at the 
credible probability of 70 %). We used the velocity range 
|Az;| < 20 km s^^ where we have X14 30 = 1-159. The 
best fit with Xmin = 0-^0 is shown in Fig. la by the sohd 
curve, whereas points with error bars give the normalized 
intensities (the corresponding S/N ~ 40). 

The best solution yields Ni = 3.9 x 10^^ cm"^ and b = 
8.8 km s~i. Of course, our analysis gives a certain range 
for allowable A^i and b. This is shown in Fig. lb, where the 
70 % confidence region for this pair of physical parameters 
is restricted by the upper part of the ellipse over the line 
& = 8 km s"'^. The second part of the ellipse under this 
horizontal dashed line may be excluded from consideration 
since the 6- value of the heaviest element Znii equals 9.3 ± 
1.3 km s""'^ (Paper I), otherwise one would have a zone 
where bn^ < 6znii- The cross in Fig. lb marks the point 
of maximum likelihood for the best fit. Two additional 
curves in Fig. la (dotted and dashed) show the L(4-0)R(1) 
profiles for the limiting values of A^i = 3.0 x 10^^ cm^^ and 
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Fig. 3. (a) - "9- function defined in Eq. vs the exci- 
tation temperature Tex- The values of 5" = 6 restrict the 
range for acceptable Tex as illustrated in panel b. (b) - 
Population density in H2 rotational levels observed in the 
Zahs = 3.3901 system. The column density Nj=i for H2 
molecules in level J — 1 per cm^ (or the upper limits of 
Nj for J = 0, 2 — 5), divided by the statistical weight, gj, 
is plotted on a logarithmic scale against the excitation en- 
ergy, Ej, in K. The corresponding values of J are shown in 
the horizontal scale at the top. The points have been fit- 
ted with one straight line using Eq.(^, and a value of Tex, 
obtained from the usual Boltzmann formula. The shown 
straight lines correspond to Tex = 50, 500, and 1400 K 



4.8 X 10 cm , respectively. Both of them have Xmin ~ 
1.06. To summarize, we set A^i = (3.9 ± 0.9) x 10" cm^^ 
and = 8.8tjj;g km s'^. 

Let us now turn to the estimation of the total A^(H2). 
As mentioned above the total molecular hydrogen column 
density depends on T^x- Typical ISM diffuse molecular 
clouds have a mean (J — 0)/(J — 1) excitation temper- 
ature of (Toi) = 77 ± 17 K, with the bulk lying in the 
range from 45 to 128 K (Savage et al. 1977). The exci- 
tation temperature Tqi reflects the relative populations 
of J = and 1 rotational levels which are controlled by 
thermal particle collisions. Therefore Tqi is approximately 
equal to the kinetic temperature, Tkin, of the cloud. The 
population of the higher rotational levels are established 
by collisions, UV pumping and radiative cascades after 
photo-absorption to the Lyman and Werner bands and 
thus Tex for the higher levels may differ from Tqi in dense 
molecular clouds (see, e.g. a review by Dalgarno 1976). 
However, for tenuous clouds, a unique temperature for all 
the levels can be found suggesting that the rotational pop- 



ulation could be dominated by collisions in a gas whose ki- 
netic temperature is near T^x- For example, the excitation 
temperature of 1120 ± 80 K toward ( Pup was measured 
for J = - 7 by Morton & Dinerstein (1976). 

There is also a tendency toward low excitation tem- 
peratures Tex ^ 100 K for the clouds with /(H2) > 0.05, 
whereas clouds with /(H2) < 0.05 tend to higher Tex 
(Spitzer & Jenkins 1975). Note that spin temperatures 
> 10^ K have been measured in DLAs with A^(Hi) > 10^^ 
cm~^ for which 21 cm line measurements are available (de 
Bruyn et al. 1996, Carilli et al. 1996). These measurements 
assume, however, a similar covering factor of the optical 
and radio background source. 

The excitation temperature can be estimated directly 
from the comparison of the population density in H2 ro- 
tational levels or indirectly from the thermal width of the 
H2 lines assuming that Tex is approximately equal to the 
kinetic temperature of the cloud. For the DLA in ques- 
tion, we know exactly that the line widths do not reflect 
the thermal motion since all species from the lightest H2 
to the heaviest Zn show approximately the same &- values 
(see Paper I). We also have not accurate measurements of 
Nj for different J. Therefore, to estimate Tex we applied 
another method. 

We choose the continuum 'windows' which are sensi- 
tive to the presence of the H2 absorption and set upper 
limits for Nj using the inequality 



, obs 

^ ^0 



(3) 



with T§^^ being the observed limit at the expected H2 
transition and the quantity calculated from (2). In 
this estimations we used &H2 = 8.8 km s"""^. The obtained 
limits are depicted by short horizontal bars in Fig. 3b. 

Now we know the upper limits iVj™ for J = 0, 2, 3, 4, 5 
and the measured A'^i value. To estimate Tex consider the 
following function : 



Vl'(Tex)=^Vj, 
J 

where for J = 0, 2, 3, 4, 5 we have 

_ M , if A^.f''^ < A^™ 
' 0, otherwise , 



(4) 



and for J = 1 

" 1, if INI^^'^ 



■07 



A^il < SN 



0, otherwise , 



where 6N = 0.9 x 10" cm^^^ 
Here Nj'^^'^ is calculated as 
B^J{J 

exp — 

91 
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Al— exp 



2-B,, 



T 



T 



(5) 



It follows from the definition of \E'(Tex) that the al- 
lowable range for T^x should correspond to the maximum 
value of 5" (here equals 6). The shape of ^E" as function 
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Fig. 4. Velocity plots of the UVES normalized data of 
QOOOO 2620 (dots and la error bars) obtained by Molaro 
et al. (2000) and the synthetic H2 spectra convolved with 
the instrumental resolution of FWHM = 6 km s~^ (solid 
curves), corresponding to the hmits with T^x = 50 K for 
J = (bottom panel) and Tex = 1400 K for J = 2-5. The 
dotted curve in the bottom panel shows the W(0-0)R(1) 
profile for Nj-i = 3.9 x 10^^ cm~^; the Doppler parameter 
equals = 8-8 km s~^ for all synthetic spectra (see also 
Fig. la) 



of Tex is shown in Fig. 3a which readily gives the range 
50 K < Tex < 1400 K. The range is rather wide since we 
have only one unblended H2 line. It should be noted, how- 
ever, that the obtained limits do not contradict the above 
mentioned observational estimations of Tex- 

Fig. 3b gives additional illustrations for our procedure. 
Here we show the measured (J = 1) and estimated (J ^ 
1) column densities against energies of the low rotational 
levels. Three different lines represent solutions for Tox = 
50, 500 and 1400 K. As seen, at low Tex the most sensitive 
are the J = and J = 1 transitions, whereas at high Tex 
the most stringent restriction is set by the populations of 
the J = 2 and J — A rotational levels. 

These two regimes are shown in Fig. 4 where dots with 
error bars are the observed normalized intensities within 
the chosen windows, and the solid curves show the H2 



synthetic profiles calculated with bn^ = 8.8 km s"-'^ and 
the limiting column densities listed in the corresponding 
panels. 

If Tex lies between 50 K and 1400 K, the total molecu- 
lar hydrogen column density A^(H2) ranges from 5.4 x 10^^ 
cm~^ to 1.7 X 10^^ cm~^ (see Fig. 5) where the lower value 
is obtained at Tox = 2i3t,. Thus we may set logA^(Il2) = 
13.98 ± 0.25. With logiV(Hi) = 21.41 ± 0.08 taken from 
Lu et al. (1996), and assuming that the molecular gas 
coincide with the Hi (i.e. filling factor is 1), one finds 
log/(H2) = -7.43 ±0.26 or /(H2) ~ 4 x lO'^. 

Since H2 molecules are formed more efficiently on dust 
grains in high Hi column density clouds, we can conclude 
that the dust number density in the Zabs = 

3.3901 ab- 
sorbing material is lower than in the ISM clouds with 
A'' (Hi) > 10^^ cm~^, provided the dust destruction rate 
is the same. 

2.3. Dust Content 

The presence of dust in DLAs is illustrated in Fig. 6 where 

the relation between log /(H2) and relative heavy element 
depletion, [Cr/Zn], is shown. In DLAs, the amount of dust 
is usually estimated from the ratio A''(Cr ll)/A/'(Znll) as- 
suming that Zn is undepletcd (e.g. Pettini et al. 1994). 

The measured data for the four extragalactic molecu- 
lar clouds (dots in Fig. 6) can be very well fitted with a 
linear law : log/(H2) ~ —5.5 [Cr/Zn] — 7.1, which means 
that these quantities are strongly correlated (dashed line 
in Fig. 6). Another four DLAs with high neutral hydro- 
gen column densities where the H2 absorption has not yet 
been detected (see Table 1) are marked by vertical arrows. 
Crosses in Fig. 6 are the Galactic interstellar medium mea- 
surements taken from the compilation of Roth & Blades 
(1995). By the horizontal lines we mark the H2 measure- 
ments in the Magellanic Clouds (Richter 2000). Unfortu- 
nately, the corresponding [Cr/Zn] data are not available 
yet, therefore, to restrict their plausible range we took the 
limiting values of [Cr/Zn] = —1.09 and —0.72 obtained 
from different sightlinc measurements in the Magellanic 
Clouds by Roth & Blades (1997). 

Qualitatively, the DLA data, including the /(H2) up- 
per limits, exhibit the same tendency as the Galactic ISM 
diffuse clouds, i.e. clouds with large amount of H2 show 
higher heavy element depletion. The linear correlation 
among the DLA systems is then a kind of upper enve- 
lope of all the different types of data points. The fact 
that the ISM measurements are found at higher Cr de- 
pletion than the DLA data is probably due to the higher 
metallicity and dust content of interstellar clouds. In fact, 
there is evidence for a correlation between dust-to-gas ra- 
tio and metallicity among DLA systems (Vladilo 1998). 
Given such a correlation, it is natural to expect that 
DLA points, which have lower metallicities than interstel- 
lar clouds, have also lower Cr depletion. The lack of any 
data points above the envelope can be readily explained 
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Fig. 5. The total molecular hydrogen column density 
A^(H2) vs Tex for a fixed value Nj^i = 3.9 x lO" cm'^ 

because clouds with low dust content (and low Cr deple- 
tion) are expected to have low molecular content. It is 
interesting that the data of the Magellanic Clouds fall in 
the range occupied by DLAs, as expected. Indeed, typical 
dust-to-gas ratios in the Magellanic Clouds are about 4 
times (LMC) and 8 times (SMC) lower as compared with 
the Milky Way (Koornneef 1982; Bouchet et al. 1985, re- 
spectively) . 

The fact that the DLA measurements lie on the en- 
velope while the upper limits lie below is puzzling. One 
possible explanation is the presence of some selection ef- 
fect that allows us to detect H2 at high Cr depletions only 
when the H2 lines are sufficiently strong. This could be 
the case if dust obscuration is important. 

The dust number density, rid, in the Zabs = 3.3901 
cloud may be estimated from the measured H2 fractional 
abundance assuming that the H2 formation rate on dust 
grains, i?, is in equilibrium with the II2 destruction caused 
by UV photons. Another channel of the II2 formation 
through the formation of the catalyst H~ (e.g., Donahue 
& ShuU 1991) may compete with the dust process in the 
partially ionized warm regions and will not be considered 
here because of the very high Hi column density in the 
2abs = 3.3901 cloud which provides most of the gas to be 
neutral. In the mostly neutral gas, H2 will be in equilib- 
rium with a density (e.g., Spitzer 1978) : 

/n(H2) = i?nn(H) , (6) 

where / is the H2 photo-dissociation rate in s"'^, n = 
n(H) -I- 2n(H2) in cm~^, and R in cm^s~^. For tenuous 
interstellar clouds, / ~ 0.11 f3o, with f3o being the photo- 
absorption rate in the Lyman- and Werner- bands. 

The H2 formation rate on dust grains can be written 
as (Spitzer 1978) : 



Rn 



(7) 



where the constant 7 ~ 0.3, the mean thermal velocity of 
the Hi atoms caused by particle collisions (wh) — 1.46 x 
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Fig. 6. Relation between H2 fractional abundance plotted 
on a logarithmic scale and relative heavy element deple- 
tion in high redshift DLAs (dots, vertical arrows label sys- 
tems from Table 1). The dashed line corresponds to the 
linear law : log/(H2) ~ — 5.5 [Cr/Zn] — 7.1 . The same 
plot for the Galactic ISM diffuse clouds based on the data 
from the compilation of Roth & Blades (1995) are shown 
by crosses. The horizontal lines mark the H2 measure- 
ments (Richter 2000) in the Magellanic Clouds (see text 
for details). References: 1,2 - Ge & Bechtold 1999; 3 - Lev- 
shakov & Varshalovich 1985, Srianand & Petitjean 1998; 
4 - this work 

lO'' V^kin in cm s~^, and {n^ad) is the mean area obscured 
by dust grains in 1 cm'^ along the sightline. 

The value of (3q depends on the local UV radiation 
field. The intergalactic background radiation at z ~ 3 
gives /3o = 2 X 10"^^ s'^ if J2i(912A) = 1 (Srianand & 
Petitjean 1998). If we take into account that locally in the 
intervening galaxy the UV emission from hot stars may 
increase Pq by a factor of k, then using (^) and (^ we 
have an estimate 



(ndfTd) ~ 5.0 X 10 



-17 



/(H2), 



(8) 



which is 



where we took n(H2)/n(H) = A^(H2)/A^(Hi) 
vahd for the QOOOO-2620 cloud. 

Given the mean Galactic value (nd(Td)iSM = 1.2 
10-"' 



^"2^ cm~i (Spitzer 1978), one has 



("-dCTd)lSM 



4.0 X 10^- 



:/(H2) 



(9) 



If K/\/rkin ^ 1 and the linear size distribution func- 
tion for the dust grains is invariant, then for the Zabs = 
3.3901 cloud the last relation yields ('T-d)DLA/("d)iSM ~ 
10-3. 
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It is reasonable to compare this quantity with the dust- 
to-gas ratio, fc, which may be estimated from any two 
elements X and Y having similar nucleosynthetic history 
(cf. Vladilo 1998) : 



Table 1. QSO candidates for H2 absorption 



k = 



10[x/H]< 



X,ISM 



/y, 



ISM 



(^10[Y/X]„. 



1 



(10) 



where the fraction in dust /x,ism refers to the Galactic 
interstellar medium. 

By applying error propagation to the column den- 
sity measurements for the FeiiAlGll, CriiA2056, and 
ZniiA2026 lines (Table 1 in Paper I), one can find that 
[Fc/ZnL= 0.03 ± 0.07 and [Cr/Zn] = 0.05 ± 0.08. Equa- 
tion (|lO| ) implies that positive values of the [Fe/Zn] and 
[Cr/Zn] ratios are unphysical because of the dust deple- 
tions /zn,isM = 0.35, /cr,iSM = 0.92, and /fc,ism = 0.94. 
Therefore we can estimate only maximum value of k tak- 
ing low bounds of [Fe/Zn] = -0.04 and [Cr/Zn] = -0.03, 
and the upper bound of [Zn/H] = —1.97. Under these con- 
ditions we have k < 1.6 x 10^'^. This value corresponds 
precisely to the dust number density ratio estimated above 
from the H2 fractional abundance. The same result implies 
that the current interpretation of the Cr/Zn ratios as due 
to differential dust depletion rather than intrinsic chem- 
ical abundance differences of Zn, as claimed by several 
authors, is correct. 

3. Future Prospects 

The accuracy of the estimated physical parameters iV(H2) 
and Tex may be improved by additional observations of 
the quasar 0000-2620 to increase the signal-to-noise ra- 
tio at the expected positions of the L(1-0)R(2) and L(9- 
0)R(4) lines where we observe rather narrow (|Aw| < 15 
km s~^) continuum windows free from the Lya forest con- 
taminations (see Fig. 4). At some level of S/N one may 
detect weak absorption features corresponding to these 
H2 transitions. For the measured ortho-H2 column density 
A^i = 3.9 X 10^^ cm~^, such detection at 3a level becomes 
possible starting from S/N ~ 60, if Tox > 10^ K (see Ta- 
ble 2). The S/N ratios listed in Table 2 were estimated in 
the following way. 

For weak absorption lines (tq ^ 1) one has (see 
LCFB) : 



W^~8.85x 10-18 iVA^/abs 



(11) 



where the equivalent width W is given in mA, N in cm^^, 
and Aq in a. 

The equivalent width detection limit, criim, defined in 
LCFB, depends on the spectral resolution, the mean S/N 
value per pixel and the accuracy of the continuum fit, 
(7c/C, over the width of the line : 



= AA 



M 



(- 



1/2 



(12) 



QSO 




7V(Hl), 


[Cr/Zn] 


r(H2) 


/"(Ha) 


0458-02 


2.04 


8.0(21) 


-0.7 


< 2.0(-6) 


5.0(-4) 


0100-^13 


2.30 


2.5(21) 


-0.24 


< 4.0(-6) 


1.6(-6) 


0112+03 


2.42 


1.0(21) 


-0.9 


< 2.0(-5) 


6.3(-3) 


1223-H7 


2.47 


3.0(21) 


-0.6 


< 9.0(-7) 


1.6(-4) 



"Upper limits from Ge & Bechtold (1999) 
''Estimated from log /(Ha) = -5.5 [Cr/Zn] - 7.1 



where M = 2.5xFWHM is the full width (in pixels) at 
the continuum level of a weak absorption line, and AA is 
the pixel size in A. 

For real data, CTc/C < (S/N)-^ since the continuum 
level is usually estimated as a mean intensity over m pixels 
within a continuum window. Thus we may write CTc/C = 
C(S/N)-i, ^^^j.^ ^ ^ l/y/f^. Then from (|l|) we have 



,1/2 



(13) 



In our calculations we set A A = 0.04 A, FWHM = 
2.5 pixels, and ^ = 0.2. After that we calculate Nj=2 
and Nj^4 using (|5|) for a given Tox and the corresponding 
W values from (^Tj). Taking = ^W, we find S/N 
values from (13). The obtained results listed in the sixth 



and seventh columns of Table 2 show that the L(1-0)R(2) 
and/or L(9-0)R(4) hnes may be detected for a reasonable 
exposure time if Tex is not very low in the Zabs = 3.3901 
absorbing cloud. 

Another issue we would like to outline in this section is 
the DLAs with ^'(Hl) > 10^^ cm-^ and Zabs > 2 from the 
compilation of Ge & Bechtold (1999). These systems are 
listed in Table 1. Similar to our case, they show high neu- 
tral hydrogen column densities but no II2 absorption in 
moderate resolution spectra (FWHM ~ 1 A). The upper 
limits for /(II2) are shown in the fifth column of Table 1. 
In the sixth column we present the expected /(II2) values 
based on the linear relationship between log/(Il2) and 
[Cr/Zn] as found in this paper. As noted in Sec. 2.3, it is 
surprising that there is only one DLA system (zabs = 2.30) 
which is consistent with our results, i.e. the predicted II2 
fractional abundance is lower than the previously obtained 
upper limit. All other systems show too low limits for 
/(H2) to be in concordance with expected values. 

One may advocate an intrinsic spread of the II2 frac- 
tional abundance among DLAs at constant [Cr/Zn] ratios 
as an explanation. It would be also necessary to investigate 
the sample of QSOs of Table 1 with data of similar reso- 
lution and S/N ratios than the one presented here. Such 
studies could help us to understand the dust formation 
history and the metal enrichment of damped QSO sys- 
tems since dust is mainly produced in the cool envelopes 
of intermediate-mass stars and in the dense shells of su- 
pernova remnants. 
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Table 2. Predicted equivalent widths for the L(1-0)R(2) 
and L{9-0)R(4) Hnes for different Tox and the S/N ratios 
which can provide 3a level detection 



2~cx 

(K) 


iVj=2 

(cm-2) 


iVj=4 

(cm-^) 


(mA) 


(mA) 


S/N 
R(2) 


S/N 
R(4) 


1400 


1.70(13) 


1.30(13) 


2.90 


7.97 


165 


60 


1000 


1.54(13) 


8.39(12) 


2.63 


5.14 


180 


90 


500 


1.09(13) 


1.80(12) 


1.86 


1.10 


260 


430 



4. Conclusions 

The outcome of our study, aiming to provide a measure- 
ment of the molecular hydrogen abundance in the low 
metallicity DLA galaxy at Zabs = 3.3901 toward QOOOO- 
2620, is as follows. 

1. In the high resolution UVES spectrum of the quasar 
0000-2620 a weak absorption feature at Aobs = 4609.4 A 
is identified with the H2 L(4-0)R(1) line from the Zabs — 
3.3901 system. 

2. The application of the standard Voigt fitting anal- 
ysis to this line gives the ortho-H2 column density of 
Nj^i = (3.9 ± 0.9) X 10" cm-2 and the Doppler width 
6h2 = 8.8tQ g km s"^ The width of the H2 line is found 
to match those for numerous unsaturated metal lines mea- 
sured at the same redshift in Paper I. 

3. The upper limits to the column densities of the W(0- 
0)R(0)-hR(l), L(1-0)R(2), L(10-0)P(3), L(9-0)R(4), L(9- 
0)R(5) lines together with the measured value of Nj^i 
lead to the excitation temperature Tex ranging from 50 K 
to 1400 K. 

4. The range of Tex and the Nj=i value yield, in turn, 
the total H2 column density of iV(H2) ~ (0.5 - 1.7) x 
10^"* cm~^ and the H2 fractional abundance of /(H2) ~ 
4 X 10^^. This is the lowest H2 fractional abundance ever 
measured. 

5. Assuming that the H2 formation rate on dust grains 
is in equilibrium with the H2 destruction caused by UV 
photons, the dust number density is estimated in the 
Zabs = 3.3901 cloud, (nd)DLA ^ 10"^ (nd)iSM- This value 
is in excellent agreement with the dust-to-gas ratio of 
k < 1.6 X 10~^ found independently from the [Cr/Zn] 
and [Fe/Zn] abundance ratios. 

6. The derived H2 fractional abundance excellently fits 
in the linear relation between log/(H2) and [Cr/Zn] mea- 
sured for other DLAs at high redshift. 

7. It is shown that future observations with higher 
signal-to-noise ratio (S/N > 60) may reveal a few addi- 
tional II2 lines from the Zabs = 3.3901 system which may 
increase the accuracy of the iV(Il2) and Tox measurements. 
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